Pause to make a column for biomass to chla ratios And also make biomass in in micro grams to more comparable to chl-a
Additive columns for combined surface and pore water nutrient as tot_ in ugL
Colored by years
Can we frame wet years as a shrinking productivity window? Or is it more a function of a mismatch in phenology ques, where peak light occurs at unsuitable streamflows (e.g. to high) for epilithic algae?
How large is that window from year to year?
How is the cumulative GPP during that time?
How is the cummulative PAR and water temp?
DO in black, metabolism in green
### big grid
WY_n_grid2 <- ggarrange(
templt,
templt_BWU,
templt_GBL,
templt_GBU,
ncol = 1, nrow = 4)
WY_n_grid2TS_plot <- ggplot(covariat_datq%>%mutate(DOY=yday(date))%>%filter(site=="BWL")%>%filter(water_year<2024), aes(x = DOY, y = PAR_surface, color =as.factor(water_year), shape = site)) +
geom_point(size = 2, alpha = 0.5, position = position_dodge(width = 0.3)) +
# geom_line(aes(x = DOY, y=PAR_inc )) +
geom_line(aes(x = DOY, y=(Q_m*100))) +
scale_color_viridis_d(option = "viridis") +
geom_point(size = 2, alpha = 0.7) +
#scale_color_manual(values = siteC_colors) +
scale_shape_manual(values = c(15, 0, 17, 2)) +
ylab(expression(PAR~at~stream~surface~(mu~mol~m^-2~s^-1))) +
theme_classic() + facet_grid(site~.)
TS_plotTS_plot1 <- ggplot(covariat_datq%>%filter(water_year<2024), aes(x = Q_m, y = GPP_mean, color =as.factor(water_year), shape = site)) +
geom_point(size = 2, alpha = 0.5, position = position_dodge(width = 0.3)) +
# geom_line(aes(x = DOY, y=PAR_inc )) +
#scale_color_viridis_d(option = "viridis") +
geom_point(size = 2, alpha = 0.7) +
#scale_color_manual(values = siteC_colors) +
scale_shape_manual(values = c(15, 0, 17, 2)) +
#ylab(expression(PAR~at~stream~surface~(mu~mol~m^-2~s^-1))) +
theme_classic() + facet_grid(site~.)
TS_plot1Max Q where GPP was obtained…
BWL : 2.763 cms
BWU : 5.990 cms
GBL : 0.622 cms
GBU : 0.299 cms
Dotted lines for onset of baseflow Dashed lines for peak PAR